Does confining the hard-sphere fluid between hard walls change its average properties? 
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We use grand canonical transition-matrix Monte Carlo and discontinuous molecular dynamics 
simulations to generate precise thermodynamic and kinetic data for the equilibrium hard-sphere fluid 
confined between smooth hard walls. These simulations show that the pronounced inhomogeneous 
structuring of the fluid normal to the confining walls, often the primary focus of density functional 
theory studies, has a negligible effect on many of its average properties over a surprisingly broad 
range of conditions. We present one consequence of this insensitivity to confinement: a simple 
analytical equation relating the average density of the confined fluid to that of the bulk fluid with 
equal activity. Nontrivial implications of confinement for average fluid properties do emerge in 
this system, but only when the fluid is both (i) dense and (ii) confined to a gap smaller than 
approximately three particle diameters. For this limited set of conditions, we find that "in-phase" 
oscillatory deviations in excess entropy and self-diflusivity (relative to the behavior of the bulk fluid 
at the same average density) occur as a function of gap size. These paired thermodynamic/kinetic 
deviations from bulk behavior appear to reflect the geometric packing frustration that arises when 
the confined space cannot naturally accommodate an integer number of particle layers. 



I. INTRODUCTION 

Confined fluids play an important role in a host of sci- 
entific phenomena and technological applications. Exam- 
ples range from the aqueous fluids that fill the cytostruc- 
tures of biological cells to the solvents that facilitate the 
operation of nano- and microfluidic devices, membranes 
for separations, and porous catalytic materials. In many 
of these systems, confinement significantly modifies the 
thermodynamic and kinetic behavior of the fluid relative 
to the bulk phase. Such modifications are generally at- 
tributed to the collective effects of the size and shape of 
the confined space and the interactions of the fluid with 
the confining surfaces. However, isolating the individual 
contributions of these various factors for study can be a 
daunting experimental task. 

Given this difficulty, one alternative approach has been 
to explore the behavior of simplified models that allow 
one to examine the implications of confinement in the 
absence of the complicating details that are present in 
experimental systems. Along these lines, a commonly 
investigated model is the equilibrium, monatomic hard- 
sphere (HS) fluid confined between smooth and parallel 
hard walls. This is arguably the most basic model that 
can capture the main entropic packing effects associated 
with fluids in confined spaces. Its characteristic inhomo- 
geneous density profile (normal to the confining walls), 
which has been a primary focus of previous investiga- 
tions, is now qualitatively understood. 1 Unfortunately, 
despite progress in elucidating some of the other prop- 
erties of this system^'2i 3 ri^£&L&£ i 10 i 1 1 a comprehensive 
picture for precisely how confinement modifies the aver- 
age thermodynamic and kinetic behavior of the equilib- 
rium HS fluid has yet to emerge. 

One of the most basic hurdles to constructing this pic- 



ture has been the lack of accurate molecular simulation 
data for the average properties of the confined HS fluid, 
a fact that may seem surprising given the apparent sim- 
plicity of the model. Ironically, the model's simplicity has 
indirectly contributed to the lack of simulation data be- 
cause it has allowed the system to be readily studied by 
approximate theories instead , 12 d 3 which are appealing 
because they are physically insightful and require only 
modest computational resources. However, the develop- 
ment of efficient algorithms for investigating systems with 
discontinuous potentials and the availability of fast com- 
puters have now made it feasible to use molecular sim- 
ulation to fully characterize the behavior of this model 
with both accuracy and precision. One aim of the present 
study is to leverage these simulation resources to take an 
important step toward completing this characterization. 

The data that we present here provide some insights 
into an important, but still poorly understood, concep- 
tual point concerning this model. Specifically, it has not 
been entirely clear how one should compare the confined 
fluid to the bulk fluid in order to elucidate the main ef- 
fects of confinement. One obvious possibility is to com- 
pare the two systems under conditions where they exhibit 
equal "average" density. The argument for choosing this 
basis of comparison is straightforward. Packing effects 
dominate the behavior of athermal systems, and average 
density is an important factor in determining how the 
particles pack. Moreover, if one controls for average den- 
sity in making the comparison, then one can hope to iso- 
late more subtle effects due to, e.g., the finite size of the 
confined system (in one direction) and the shape of the 
density profile (i.e., the "layering"). Another possibility 
is to compare the two systems at equal activity, where 
the bulk and pore fluids exhibit different average densi- 
ties. The advantage in doing so is also obvious. Equality 
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of activity is a relevant experimental constraint on the 
chemical equilibrium that is established between the bulk 
and pore fluids. 

The complication in comparing the two systems at the 
same density is that there are two different definitions for 
average density that are commonly invoked: p — N<r 3 /V 
and ph = Na 3 /Vh- Here, TV refers to the number of par- 
ticles, and a is the particle diameter. The difference be- 
tween the two is that V — AH is the total volume of the 
confined fluid (i.e., A is the area of a wall in contact with 
the fluid, and H is the distance from "wall surface to wall 
surface" ) , while Vh — Ah is the smaller volume accessible 
to the particle centers (i.e., h = H — a). While densi- 
ties based on these two definitions converge in the limit 
H — > oo, they can be quite different for severely confined 
fluids. We are not aware of any systematic comparisons 
for how the thermodynamic properties of this system de- 
pend on p and ph, respectively. 

However, even in the absence of such comparitive stud- 
ies, it is easy to imagine that one might indeed arrive at 
qualitatively different conclusions about the implications 
of confinement depending on whether p or ph is chosen 
as the basis for comparison. To appreciate this point, 
consider that ph diverges in the limit where the gap size 
H is reduced, at fixed N/A, to the size of one particle di- 
ameter a (i.e., the two-dimensional fluid limit), whereas 
p and many other fluid properties of interest remain fi- 
nite. This type of consideration alone hints that p might 
be the more suitable density variable of the two for mak- 
ing comparisons to the bulk fluid, and indeed p naturally 
emerges in the thermodynamic analysis of confined HS 
fluids^ 

More concrete evidence supporting the use of p 
rather than ph for comparing confined and bulk flu- 
ids comes from studies of transport properties. Specif- 
ically, it has recently been demonstrated via molecular 
simulation^ that the self-diffusivity of the confined HS 
fluid parallel to the confining walls, over a broad range 
of equilibrium conditions, is very similar to the diffusion 
coefficient of the bulk HS fluid if the two systems are 
compared at the same value of p. In other words, the 
specific details of the inhomogeneous packing structures 
have only minor influence on the average single-particle 
dynamics of the confined fluid, as long as one controls for 
the average overall density p. Alternatively, if one instead 
compares the behaviors of the bulk and confined systems 
at equal values of ph, one arrives at the conclusion that 
confining a HS fluid between hard walls has the effect of 
significantly speeding up its dynamics. This latter ar- 
tificial conclusion is related to the fact that N/A must 
vanish if ph is to remain constant in the limit H — ► a. 
As a result, even if the numerical value of pt is chosen to 
be indicative of a dense bulk fluid, the actual average in- 
terparticle separation and particle mobility in the lateral 
direction will generally be very large (e.g., comparable to 
a dilute gas) when the fluid is confined to small enough 
H. 

In this paper, we follow up on some of these initial ob- 



servations by presenting a more comprehensive study for 
how fluid density and confinement (between hard walls) 
affect the thermodynamic and kinetic properties of the 
HS fluid. We broadly focus our investigation on four 
main questions. The first pertains to the equation of 
state of the confined fluid (i.e., how the average trans- 
verse and normal components of its pressure tensor vary 
with average density). Specifically, we are interested in 
how the behaviors of these pressure components depend 
on the volume definition invoked, i.e., V versus Vh- Does 
use of either defintion produce relationships similiar to 
the equation of state the bulk HS fluid? Second, what 
are the effects of confinement and average density on the 
transverse self-diffusivity of fluids confined to pores nar- 
rower than those previously examined^ (i.e., H < 3.5er)? 
Third, how do the behaviors of the confined and bulk 
HS fluid systems compare under conditions of equal ac- 
tivity as opposed to equal density? Finally, does the 
robust relationship between excess entropy s ex (relative 
to ideal gas) and self-diffusivity D, previously discovered 
for fluids confined to gap sizes larger than H = 3.5er 
in this system^ continue to hold for very narrow pores 
(H < 3.5er)? By addressing these four questions, we can 
make significant headway not only in indentifying the re- 
gions in the H — p and H — £ planes of parameter space 
where the confined HS system significantly deviates from 
the bulk HS fluid, but also in probing the microscopic 
mechanisms for such deviations. 



II. SIMULATION METHODS 

To explore these issues, we have calculated the ther- 
modynamic properties of confined and bulk HS fluids 
using grand canonical transition-matrix Monte Carlo 
(GC-TMMC) simulations ^±2. and we have tracked their 
single-particle dynamics via discontinuous molecular dy- 
namics (DMD) simulations^ To simplify the notation 
in this article, we have implicitly non-dimensionalized all 
quantities by appropriate combinations of a characteris- 
tic length scale (which we take to be the HS particle di- 
ameter a) and time scale (which we choose to be cry/m{3, 
where m is particle mass, (3 = [k-^T]^ 1 , fee is the Boltz- 
mann constant, and T is temperature). As a result, all 
quantities with dimensions of energy are understood to 
be "per fceT", the only energy scale in the problem. 

The DMD simulations each involved N — 1500 identi- 
cal HS particles. For the bulk fluid, the particle centers 
were contained within a cubic simulation cell of Vh = 
N/ ph, and periodic boundary conditions were applied in 
all three directions. For the confined fluid, particle cen- 
ters were contained within a rectangular parallelepiped 
simulation cell of Vh = h x h y h z , where h z = H — 1 and 
h x = hy = [N/ (hzph)] 1 ' 2 ■ Periodic boundary condi- 
tions were applied in the x and y directions and per- 
fectly reflecting, smooth hard walls were placed so that 
particle centers were trapped in the region < z < h z . 
The self-diffusivity D of the fluid was obtained by fit- 
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ting the long-time (t > 1) behavior of the average mean- 
squared displacement of the particles to the Einstein rela- 
tion (Ar 2 ) = 2dDt, where Ar 2 , corresponds to the mean- 
square displacement per particle in the d periodic direc- 
tions (d = 2,3 for the confined and bulk fluid, respec- 
tively). To verify that system-size effects in the periodic 
directions on D were insignificant, we checked that our 
calculated values for D for several state points compared 
favorably with those we obtained using either N = 3000 
or N = 4500 particles. 

The GC-TMMC simulations each utilized a simulation 
cell of size Vh = 1000. For the bulk fluid, the cell was cu- 
bic with h x = h y = h z = 10. For the confined fluid, the 
cell was a rectangular parallelepiped with h z — H — 1 and 
h y = h z = ^1000//i 2 . GC-TMMC simulations require 
a specified value for the activity §il (i.e., N is allowed 
to fluctuate), and we set f = 1 in all cases. The key 
quantities that we extracted from the simulations were 
the normalized total particle number probability distri- 
bution II(iV) and the iV-specific spatial density distribu- 
tion p(N, r) , both evaluated over a range of particle num- 
bers spanning from N — to N = 984. Thermodynamic 
properties at other values of activity £ were readily ob- 
tained via the histogram reweighting technique 18 to shift 
the original H(N) distribution to one representative of 
the particle numbers visited at the selected £. We found 
that we obtained statistically indistinguishable results for 
systems with Vh = 500, indicating again that noticeable 
artifacts associated with system size were not present. 

By employing basic arguments from statistical 
mechanics one can use the equilibrium informa- 

tion from GC-TMMC simulations to compute thermo- 
dynamic properties of interest. Specifically, the grand 
potential f2 can be calculated directly from the normal- 
ized particle number distribution ; 19 ! 20 

Q = inn(o). (l) 

For the bulk HS fluid, we also have V — Vh, and thus 
p = J^n NIl(N)/V — ph- Moreover, the pressure of the 
bulk fluid P is equal to the negative of the grand potential 
density, P — —Ct/V. On the other hand, for the HS fluid 
confined between hard walls, we have V = Vh/(1 — H^ 1 ), 
and thus p = J2n nu ( n )/ v = C 1 - H ~ X )Ph- In this 
case, negative grand potential density — f2/V represents 
an average transverse pressure acting parallel to the con- 
fining walls. 21 In the reduced units adopted here, the 
component of the pressure tensor acting normal to the 
walls is equal to the local fluid density in contact with a 
hard wall, P Z (N) = p(N,z = 0.5) = p(N, z = H - 0.5), a 
consequence of an exact statistical mechanical sum rule 
for this system^ Finally, the molar excess entropy s cx = 
5 ex /iV is determined using the following expression; 7 ^ 

S CX {N) = ln[n(iV)/n(0)] - JV1< + In AH 

-N\nN + J p(N, r) In p(N,r)dr. (2) 

Below, we describe how the above methods were em- 
ployed in this study to characterize the behaviors of the 



confined and bulk HS fluids. 



III. RESULTS AND DISCUSSION 
A. Volume definition and the equation of state 

One of the most practically important and well- 
understood properties of the bulk, equilibrium HS fluid is 
its equation of state P{p), which quantifies how its pres- 
sure varies with density. For densities below the freez- 
ing transition (p < 0.943), this relationship is accurately 
described by the semi-empirical Carnahan- Starling equa- 
tion P(p) « Pcs(p) = p(l + ^+^ 2 -0 3 )/(l-0) 3 r^ where 
4> = np/6 is the packing fraction of the spheres. 

Much less is known about the global behavior of the 
pressure tensor for the HS fluid confined between smooth 
hard walls. One obvious question is, do the relationships 
between the transverse and normal components of the 
pressure tensor and "average" density (defined as either 
p or ph) show quantitative similarities to the equation 
of state of the bulk fluid? Although the inhomogeneous 
structuring of the fluid might be expected to give rise to 
some nontrivial deviations from bulk fluid behavior, the 
main qualitative trends should be the same: compress- 
ing the fluid increases the interparticle collision rate and, 
consequently, the individual components of the pressure 
tensor. 

In Fig. [H we compare the bulk fluid equation of state 
to our GC-TMMC simulation data for the average trans- 
verse and normal components of the pressure tensor. We 
focus here on confined fluids with H — 3.5, 6, 8.5, and 
16. In top panels (a) and (b), negative grand potential 
density is plotted versus average density, adopting the 
Vh (center accessible) and V (total) volume conventions, 
respectively. The density dependencies of P z are simi- 
larly displayed in panels (c) and (d). Focusing on plots 

(a) and (c), one finds a family of curves for —Cl/Vh and 
P z that are qualitatively similar to the bulk fluid be- 
havior, with the main difference being that systems with 
smaller H have weaker ph dependencies (higher apparent 
compressibilities). This difference appears logically con- 
sistent with the earlier observation 7 that confined HS 
fluids also have faster single-particle dynamics as com- 
pared to the corresponding bulk fluid with the same ph- 

Interestingly, the corresponding quantities plotted in 

(b) and (d) using the total volume V convention ap- 
proximately collapse onto a single curve. This means 
that the p dependencies of both — Q,/V and P z , for each 
value of H investigated, can approximately be described 
by the equation of state of the bulk fluid P(p). This 
trend also appears consistent with the approximate col- 
lapse of self-diffusivities for confined HS fluids onto the 
bulk behavior when plotted together on a single graph 
versus p.- Although there are clearly some quantitative 
deviations from bulk behavior for the smallest pores in 
panels (b) and (d) of Fig. [TJ we found that the follow- 
ing simple relationship can describe the p dependence of 
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FIG. 1: Equation of state of a confined HS fluid between hard 
walls separated by center-accessible distance H — 1. Top pan- 
els show the negative of the grand potential density (average 
transverse pressure) versus average fluid density calculated 
using (a) the center-accessible volume Vh and (b) the total 
volume V. Bottom panels illustrate the normal pressure ver- 
sus average density calculated using (c) the center-accessible 
volume Vh and (d) the total volume V . 



the grand potential density to within at least 25% for 
H > 3.5: 



n( P ,H)/v^-p cs (p). 



(3) 



We will use this approximate relationship below to help 
construct an analytical model for predicting the excess 
adsorption of fluid in a model slit pore. 



B. Interfacial free energy and excess adsorption 

Given the approximate collapse of the thermodynamic 
data for the confined HS fluid when plotted against av- 
erage density p, it is natural to ask whether there is a 
connection to the behavior of the interfacial free energy 
and the surface excess adsorption of the fluid at a single 
hard wall. 



The interfacial free energy of the HS fluid near a hard 
wall is defined as the excess grand potential of the fluid 
(relative to bulk) per unit fluid-wall contact area. Similar 
to average density, its numerical value depends on the ar- 
bitrary choice of dividing surface ) 25 ' 26 although different 
choices provide equivalent thermodynamic descriptions 
of the system if applied self-consistently. If one chooses 
the plane of closest approach of the particle centers to the 
wall as the dividing surface, then the following expression 
yields the interfacial free energy: 



7,f = lim 



Vh 



(H-l) 



(4) 



where Pb is the pressure of the bulk fluid in equilibrium 
with the pore fluid. Stated differently, p of the pore fluid 
is determined by H and the requirement that it adopt the 
same activity £ as the bulk HS fluid of pressure Pb. There 
is an accurate approximate equation due to Henderson 
and Plischke 2 ^ for predicting how 7^ depends on the 
packing fraction of the bulk fluid 0b = np^/d, 
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[1 + (44/35)0 b - (4/5)0g 
(Wb) 3 



(5) 



If one instead chooses the physical surface of the wall to 
be the dividing surface, then a slightly different equation 
emerges: 
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lim 

fi->oc 



V 



Pb 



H 
~2 



lh + ^b/2 



(6) 
(7) 



Substituting the Carnahan-Starling equation of state for 
Pb and Eq. [5] for 7^° into Eq. [7] results in the following 
analytical estimate for 7 00 , 



7 



3 , [l-(l/2)0 b -(31/35)0 2 + (l/5)0 b ] 
'b — ( 8 ) 
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(1 - M 3 



Given that we have already observed that other proper- 
ties of the confined HS fluid approximately collapse when 
plotted versus p (based on total volume V) , we choose to 
focus our attention from this point forward on 7 00 , the 
interfacial free energy that is also based on V. 

As we demonstrated in the previous section, one can 
readily determine the quantities on the right-hand side 
of Eq. E] for hnite values of H using GC-TMMC simula- 
tions. As a result, these simulations might also provide a 
reasonably accurate means for estimating 7 00 , assuming 
that H can be chosen large enough so that the pertur- 
bations to the fluid caused by the two confining walls do 
not significantly interfere with one another (i.e., so that 
so-called "finite-size" or frustration effects of confinement 
do not occur). Although, it is not clear a priori how large 
H must be to achieve this, one might reasonably expect 
that the pore would need to be at least several particle 
diameters in width. 

As a test of this idea, we present in Fig. f2] values 
of the quantity [fl(p, H)/V + P b ] (H/2) calculated from 
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FIG. 2: The quantity (H/2) [Q/V + P b ] calculated from our 
GC-TMMC simulations for various H along with the H — > 
oo limit, 7 00 , computed using Eq. [5] Data are plotted as a 
function of 0b, the packing fraction of the bulk HS fluid that 
is in equilibrium with the pore fluid. 




our GC-TMMC simulations for various H along with the 
single-wall quantity 7 00 of Eq. [5J which is the H — > 00 
limit. All data are plotted as a function of bulk pack- 
ing fraction 0b. Interestingly, the plot reveals that the 
simulated curves for H > 3.5 all collapse, to within an ex- 
cellent approximation, onto that for 7 00 . In other words, 



r 



n(p, h) 
v 



H 
~2 



(9) 



independent of H for H > 3.5. This implies that single- 
wall behavior such as 7 00 can be estimated with great 
accuracy in this system from a simulated slit-pore of 
width H = 3.5, which can only accomodate a fluid film 
three particle layers thick. This result, while very robust, 
is somewhat surprising because the single-wall density 
profiles decay slowly enough to expect appreciable in- 
terference or frustration effects at pore sizes as small as 
H = 3.5. However, similar to the picture that emerged 
from the behavior of the equation of state in the previous 
section, any interference that does occur apparently can- 
cels in determining the average properties of the confined 
HS fluid, which remain remarkably "bulk-like" even for 
these very thin films. 

So, when do interference effects due to packing frus- 
tration of the wall-induced particle layers begin to oc- 
cur? We can probe this issue by taking the analysis one 
step further. Specifically, if one uses Eq. [3] to substitute 
for H)/V in Eq. [3J differentiates both sides of Eq. M 
with repect to chemical potential, and invokes the Gibbs 
adsorption equation dj°°/dp = — r°°, then upon rear- 
ranging one arrives at the following simple equation for 
predicting the pore density p: 



P~ Pb + 



2T C 



H 



(10) 



The quantity T°° is the standard surface excess density 
for a HS fluid next a to single hard wall, and, within the 



FIG. 3: Pore fluid density p as a function of pore width H 
at different values of activity [ln£ = 0.4, 4.4, and 8.8]. Filled 
symbols and solid lines correspond to the GC-TMMC data 
and the predictions of Eq. 1101 respectively. The dashed lines 
correspond to the bulk density ph for a given activity £. 



above approximations, it is given by 

oo _ 30b [l + ai0b + Q20b + a 3(f>b + a 40b] 

" tt(i + + m - 403 + 4>i) 



(11) 



where a\ = 1, a-2 = —221/70, 03 = 4/5, and 04 = —1/5. 

Fig. [3] shows the predictions of the simple analytical 
model of Eq. \W\ and [IT] compared to the simulated pore 
density p as a function of H. ^From the plot, it is 
evident that the average pore density can be predicted 
based on knowledge of only the single- wall surface excess 
r°° unless the fluid is both dense and confined to pores 
narrower than approximately three particle diameters. 
Under those restrictive conditions, the single-wall model 
misses the emergence of oscillations in the pore density. 
These oscillations cannot be solely attributed to single- 
wall "layering" in the density profile because pronounced 
layering also occurs for dense fluids with H 3> 3, where 
the analytical model is still very accurate. Rather, the 
oscillations must be due to packing frustration associated 
with the interference of the layers emerging from the two 
confining walls, which apparently becomes significant in 
this system only for very narrow pores and high fluid 
density. 



C. Comparing bulk and pore fluid self-diffusivities 

The last two sections demonstrated that some of the 
average thermodynamic properties of the confined HS 
fluid are very similar to those of the bulk fluid if the 
two systems are compared at equal values of the average 
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FIG. 4: Self-diffusivity D as a function of pore width H at 
different pore fluid packing fractions (f>. For <f) = 0.45, crosses 
mark regions for which the confined system penetrates into 
the fluid-solid coexistence region or the solid phase region on 
its equilibrium phase diagram.— 



density p (based on the total system volume). Deviations 
occur only when the fluid is both dense and confined to 
pores narrower than approximately three particle diam- 
eters. In a previous study— we have also shown that 
the self-diffusivity D of the confined HS fluid is approx- 
imately equal to that of the bulk fluid with the same p 
for H > 3.5 over a fairly broad range of p. Here, we care- 
fully investigate the _ff-dependency of pore self-diffusivity 
at constant p for narrow pores, with a focus on under- 
standing when packing frustration causes the correlation 
between D and p to break down. We also investigate the 
-H-dependency of D for the confined fluid under the con- 
straint of constant imposed activity £. We find that this 
latter behavior can be essentially predicted in advance, 
given the known connection between D and pL and the 
ability to predict p from £ and H discussed in the previ- 
ous section. 

We begin here by examining how H affects D at con- 
stant p using the DMD simulations described earlier. 
Specifically, we plot in Fig. [4] the self-diffusivity D of 
the bulk and confined HS fluid for H — 2 to 5 and vari- 
ous pore packing fractions (<fr = irp/6 = 0.15, 0.30, 0.40, 
and 0.45). What is plainly evident is that up to fairly 
dense packing fractions (<fi < 0.40), D of the confined 
fluid shows no significant deviations from bulk behavior 
(dashed line) even when in very restrictive pores (e.g. 
H = 2). In fact, quantitative deviations are prominent 
(> 25%) only in the high density (<fr > 0.4) and small pore 
(H < 3) limit. Note that an equilibrium fluid at <j) — 0-45 
cannot be accessed over the full H range because the sys- 
tem penetrates into the fluid-solid coexistence region or 
the solid phase region on its phase diagram^ 

To gain a more physical understanding of the varia- 
tions in D at constant <f> that occur under conditions of 
high (f> and low H, we plot in Fig.[5]the 2D projections of 
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FIG. 5: 2D projections of typical instantaneous particle con- 
figurations of the confined HS fluid are shown (top) along with 
equilibrium density profiles p(z) (bottom), where z represents 
the positional coordinate normal to the walls. 



instantaneous particle configurations of the confined HS 
fluid for H = 2.0, 2.4, and 3.0 at <f> = 0.40, state points 
that show very different dynamical behaviors. We also 
present the corresponding density profiles p(z) normal 
to the walls. This figure shows well-developed layering 
structures for both H = 2.0 (two particle layers) and 
H = 3.0 (three particle layers). However, the system at 
H = 2.4 shows considerably more packing frustration. 
In particular, the individual density peaks are reduced 
in this case because the spacing is such that it is "in be- 
tween" distances that naturally accomodate either two 
or three layers. The pore diffusivity is also lowest for 
4> = 0.4 at H = 2.4 as shown in Fig. [4] Similar os- 
cillations in D, which are much smaller in magnitude 
and decay with increasing H, occur at larger separations 
with the minima again coinciding with spacings that do 
not naturally accommodate an integer number of par- 
ticle layers. In short, for small enough pores and high 
enough densities, the frustrated layering of particles nor- 
mal to the confining walls significantly slows down the 
single-particle dynamics in the direction parallel to the 
walls. 

The trend that increased layering leads to faster dy- 
namics may initially appear counter-intuitive, especially 
if one tries to understand it by drawing an analogy with 
the bulk HS system. In the bulk HS system, compress- 
ing the fluid increases the structural orderin g 28 ' 29 but 
reduces the self-diffusivity. In contrast, as can be clearly 
seen in Fig. 3] and [5l increased layering in the normal di- 
rection (i.e., less uniform density profiles) correlates with 
faster dynamics. However, these two represent funda- 
mentally different systems undergoing different changes. 
In the bulk HS system, increasing the density not only 
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increases the structural order, but it also reduces the en- 
tropy (or average free volume) of the particles in the fluid. 
This compression-induced reduction in free volume is not 
surprisingly correlated with slower dynamics. However, 
the confined HS system actually maximizes its entropy 
(or average free volume) at fixed average density and H 
by adopting an inhomogeneous density profile with pro- 
nounced layering^ Our results show that, for constant 
p, the values of (small) H that frustrate the ability of 
the system to form an integral number of particle lay- 
ers also tend to reduce the single-particle mobility in the 
direction parallel to the walls. We return to investigate 
the potential connection between dynamics and entropy 
of the confined HS fluid in the next section. 

Another important point concerning the frustration- 
induced oscillations in D of Fig.0]is that they are distinct 
from the oscillations in D that occur as a function of 
H at constant activity £ i 12 > 13 The latter are inevitably 
impacted by oscillations in average pore density, whereas 
the average density is being controlled for (held constant) 
in Fig. [H Deviations from bulk behavior at fixed average 
density are purely frustration-induced finite-size effects, 
and the relative importance of these types of deviations 
has been a long-standing question in the study of confined 
fluids^- 

Interestingly, if one compares the locations of the oscil- 
lations in D versus H at <j) = 0.40 in Fig. Q] with the fluid- 
solid phase boundary of this system presented by Fortini 
et al.£ one also finds a strong correlation between slow 
dynamics and proximity of the fluid to the phase bound- 
ary. In other words, the same packing frustration that is 
giving rise to slow dynamics also appears to ultimately 
promote the formation of an ordered solid phase. This 
argues that the effect of confinement on the phase dia- 
gram of the system can provide important insights into 
how confinement impacts single-particle dynamics. The 
consequences of this could be significant for the strate- 
gies that are typically employed to study supercooled 
and confined liquids. For example, weak polydispersity is 
commonly incorporated into model fluid systems in order 
to study them under conditions where the corresponding 
monatomic fluid would rapidly crystallize. A cautionary 
note that follows from the above discussion is that one 
should not readily assume that the behaviors of the poly- 
disperse and monatomic systems are trivially related, and 
that the former only differs from the latter in that its liq- 
uid state is kinetically accessible over a broader range of 
conditions. The phase diagrams of polydisperse materials 
are considerably more complex than monatomic systems 
(even in the bulki^) Thus, one should expect confine- 
ment to impact the dynamics of polydisperse systems in 
ways that are not easily relatable to the behavior of the 
corresponding monatomic fluids. 

We now turn our attention to the independent diffu- 
sivity behavior of the confined HS fluid at fixed activity 
£ (i.e., in chemical equilibrium with the bulk). As can 
be ascertained from the strong correspondence between 
D and p in Fig. [?J the dynamical behavior at constant 




H 

FIG. 6: Average density p and self-diffusivity D as a function 
of pore size H for the confined HS fluid in equilibrium with the 
bulk HS fluid at at a given activity £ [ln£ = 0.4, 4.4, and 8.8]. 
Dashed lines correspond to the density pb and self-diffusivity 
D of the bulk HS fluid at the given activity £. 

£ can be largely predicted in advance if one simply has 
knowledge of how H influences p at constant £ (e.g., from 
simulation or the analytical model of Eq. [TO] and ITT]) . In 
Fig. we provide the independent data along constant 
£ paths for the quantities p and D determined from GC- 
TMMC and DMD simulations, respectively. One initial 
observation is that p is always less than pi, for finite H, 
and, as should be expected based on this, D is larger 
in the pores than in the equilibrium bulk fluid. Note 
that this type of physically-intuitive connection between 
average density and dynamics would be completely lost, 
however, if one instead chooses ph as the definition for 
average density, which is significantly greater than ph for 
finite H . More generally, the reliability of approximate 
theories for transport properties in inhomogeneous fluids 
could be particularly sensitive to how averaging is han- 
dled, which might help to explain why an earlier kinetic 
theory — predicts that confining a fluid at constant £ sig- 
nificantly decreases D, the opposite of what is seen in the 
MD simulation data of Fig. [6l 

A second observation about the data in Fig. [6] is that 
there are negative oscillatory deviations in p (relative 
to bulk) with H at high £ in small pores, which one 
might expect to produce similar positive oscillations in 
D. However, the frustration-induced negative deviations 
from bulk behavior in the D versus H relationship at 
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FIG. 7: (a) Self-diffusivity D and (b) excess entropy s ex as a 
function of pore size H for the confined HS fluid at a given 
pore packing fraction (f> — ivp/6. Data from top to bottom 
correspond to (f> = 0.15, 0.3, and 0.4. (c) Self-diffusivity D 
and (d) excess entropy s ex as a function of pore size H for 
the confined HS fluid at a given activity f. Data from top to 
bottom correspond to ln£ = 0.4, 4.4, and 8.8. 



constant p shown in Fig. [5] appear to largely cancel this 
effect. The net result is that D is strikingly similar to 
bulk behavior, even for small H, along paths of constant 
(and sufficiently high) £. 



D. Diffusivity and excess entropy 

The oscillatory data in Fig. Q] clearly show that aver- 
age density alone cannot predict the self-diffusivity of the 
HS fluid if the fluid is both dense and confined to a pore 
smaller than approximately three diameters. Is there an- 
other thermodynamic quantity that can predict diffusiv- 
ity behavior is these narrow pores? One promising can- 
didate is the excess entropy s CK (relative to ideal gas), 
which recent DMD simulations? demostrate, to an excel- 
lent approximation, determines the self-diffusivity of the 
HS fluid confined between hard walls for H > 3.5. Here, 
we explore its relationship to self-diffusivity in smaller 
pores. 

Fig [7] shows the data for D and s cx of the confined 
fluid collected from our DMD and GC-TMMC simula- 
tions [data at fixed pore packing fraction ip provided in 




FIG. 8: Self-diffusivity D vs negative excess entropy per par- 
ticle -s ex for the bulk HS fluid (solid curve) and for the HS 
fluid confined between smooth hard walls (symbols). The 
filled symbols correspond to the confined system at a fixed 
pore packing fraction cj> as shown on the legend and empty 
symbols are for system at a given activity £ for which bulk 
packing fraction is given on the legend. 



panels (a) and (c), and data at fixed activity £ provided 
in (b) and (d)] . Irrespective of the thermodynamic path, 
strong qualitative correspondence is observed between D 
and s cx , including the "in-phase" oscillations that emerge 
at small H. In other words, self-diffusivity and excess 
entropy appear to be affected in a very similar way by 
confinement, even for the very narrow pores. 

To scrutinize the quantitative accuracy of the relation 
between the two variables, we also plot all data corre- 
sponding to constant pore packing fraction <\> (filled sym- 
bols) and constant £ (empty symbols) paths in Fig [8] in 
the D-s cx plane. As can be seen, most of the data falls 
very close the curve for the bulk HS fluid, indicating 
that excess entropy (a static quantity) can indeed ap- 
proximately predict the implications of confinement for 
self-diffusivity. The largest deviations are for the fluid 
that has the highest pore packing fraction of <f> = 0.4. 

This data is yet one more manifestation of a larger 
trend seen throughout this paper. Namely, that the con- 
fined HS fluid, by measure of many of its average prop- 
erties, has behavior very similar to that of the bulk fluid. 
It changes character only under a fairly restrictive set of 
conditions, when the pore fluid is dense (<j> > 0.4) and 
when it is confined to pores smaller than approximately 
three particle diameters in width. 
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IV. CONCLUSIONS 
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In conclusion, we have presented new comprehensive 
simulation results for the HS fluid confined between 
smooth hard walls. The results elucidate thermodynamic 
and dynamic behavior of this system over a wide range 
of system conditions, further clarifying the precise role 
of confinement on average fluid properties and the most 
useful way to define average density for this system. One 
perhaps unexpected result is that, for most conditions, 
the average behavior of the confined HS fluid is very 
similar to that the bulk fluid. Frustration-induced fi- 
nite effects do emerge in this system, but they are only 
prominent for very small pores (dimensions smaller than 
approximately three particle diameters) and high fluid 
densities where the system approaches the confinement- 
shifted fluid-solid phase boundary. 
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